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ABSTRACT 

We study the effect of the relative velocity of dark matter and baryonic fluids after 
the epoch of recombination on the evolution of the first bound objects in the early 
universe. Recent work has shown that, although relative motion of the two fluids is 
formally a second order effect in density, it has a dramatic impact on the formation and 
distribution of the first cosmic structures. Focusing on the gas content, we analyze the 
effect of relative velocity on the properties of halos over a wide range of halo masses 
and redshifts. We calculate accurately the linear evolution of the baryon and dark 
matter fluctuations, and quantify the resulting effect on halos based on an analytical 
formalism that has been carefully checked with simulations in the case with no relative 
velocity. We estimate the effect on the abundance of and gas fraction in early halos. 
We find that the relative velocity effect causes several changes: (i) the characteristic 
mass that divides gas-rich and gas-poor halos is increased by roughly an order of 
magnitude, from 2 x Mq to about 2 x 10^ M©; (ii) this characteristic mass has 
a large scatter (full width at half maximum is ~ 1.5 x IO^Mq at z = 20); (iii) the 
fraction of baryons in star-less gas minihalos is suppressed by a factor of 4 at z = 20; 
(iv) the fraction of baryons in halos that can cool and form stars is suppressed by a 
factor of 1.5 at z = 20; and (v) there are enhanced spatial variations of these various 
fractions. 

Key words: dark ages, reionization, first stars, large scale structure of the Universe. 
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One of the most important questions in astrophysics today is 
understanding the formation and evolution of the first bound 
structures. Significant theoretical and observational efforts 
are devoted to understanding the properties of the first 
galajcies and minihalos, at what redshifts they form and how 
they influence the epoch of reionization. Observations, most 
notably of the cosmic microwave background (CMB), have 
established the basic parameters for the in itial conditions for 
structure formation (jSennett et afllioOSl l. thus providing a 
foundation for theoretical work on the first structures. Ad- 
vances in computation have made it possible to simulate the 
formation of the f i rst stars |Abel. Brvan fc NormanI l2002l : 
lO'Shea et afllixil : lYoshida et al.lbOOs !'). Meanwhile, several 
efforts are underway to probe the structure of the intergalac- 
tic medium (IGM) during the reionization epoch using the 
21 cm line of hydrogen, and second-generation experiments 
may be able to explore the early stages of reionization. 
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To answer these and many other questions it is imper- 
ative to know the correct initial conditions that led to the 
formation of the first bound objects and to account for all 
subtle effects that influence evolution of the density per- 
turbations after recombination. The critical role of initial 
conditions has been discussed bv lNaoz. Yoshida. fc Barkanal 
(2010), who showed that three commonly used setups lead to 
significantly different abundances and properties of the first 
star-forming gas clouds as well as first gas-rich minihalos. 

There are two major classes of early-type objects that 
must be analyzed. The first class consists of large halos in 
which the gas can cool and form stars; these are the pre- 
sumed sites of the first dwarf galaxies, which represent the 
first source of metals in the Universe, and provide ultra- 
violet photons that begin the dec oupling of the hydrogen 
spin temperature from the CMB (|Madau. Meiksin fc Reea 
'1997*) and eventually start the epoch of reionization. The 
second class consists of smaller halos ("minihalos") that 
are too small for molecular cooling, but still affect the 
epoc h of reionization by acting as a sink for ionizing pho- 
tons (jHaiman. Abel fc Madaull200ll : iBarkana fc Loebll2002l : 
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llliev. Scannapieco fc Shapiro! l2005l : ICiardi et al.l 120051 ) and 

may gen erate a 21 cm signal from collisional ex citation of 
Hi (e.g. llliev et~aLl l2003l : TFurlanetto fc OhlbOOel ). It is im- 
portant to understand both the abundance and distribution 
of halos, as well as the precise boundaries separating halos 
that undergo cooling and star formation, those that collect 
baryons in their potential wells but do not cool, and the 
lightest halos that exist only as dark matter structures and 
do not collect gas. 

An important effect that was previously overlooked is 
tha t of the relative velocity of dar k matter and baryonic flu- 
ids (|Tseliakhovich fc Hiratal2010l ). This effect leads to power 
suppression on scales that correspond to the first bound ha- 
los between 10* Mq and 10* Mq and delays the formation 
of the first objects. More importantly this effect introduces 
scale-dependent bias and stochasticity, leading to significant 
qualitative changes in the distribution of the first objects. 
The relative velocity effect is especially important on the 
small scales where the first stars and galaxies form. Intro- 
duction of this effect dramatically changes the gas distri- 
bution inside the first halos and changes the char acteristic 
mass of gas-rich objects. lOalal. Pen fc SeliakI (|2010l ) recently 
calculated analytically the effect on the gas content of halos 
and found a large effect on the fluctuations of the Lyman-a 
background at high redshifts. Their analysis, however, was 
based on a very simplified model of which halos can form 
stars and in what abundance. In this paper we carry out a 
detailed analytical study of the distribution of gas and stars 
in the first halos. 

The rest of the paper is organized as follows. Section [2] 
reviews the re lative velocity effect (Sec. 12.111 and improves 
the analysis of lTseliakhovich fc Hirata l|2010! ) to account for 
spatial variation of the sound speed fSec. l2.2) l. Section[3]then 
investigates the early halos and their gas content, focusing 
on computation of the filtering mass (Sec. 13. ip and then ex- 
amining the fraction of baryons in minihalos and in larger 
halos that can cool, including an analysis of spatial varia- 
tions in the baryon budget (Sec. l3.2l and 13. Sp . We summarize 
our results in Sec.|4]and compare them to other recent work. 

The numerical results and plots shown in this pa- 
per assume a cosmology with present-day baryon density 
r2b,o = 0.044, CDM density = 0.226, dark energy den- 
sity r^A.o = 0.73, Hubble constant Ho = 71 kms"^Mpc"\ 
and adiabatic primordial perturbations of variance A^(fc^) = 
2.42 X 10"^ at it* = 0.002 Mpc"^ and with slope Us = 0.96. 



2 INITIAL CONDITIONS FOR HALO 
FORMATION 

In this section we detail the formalism necessary for gen- 
eration of correct initial conditions taking into account two 
important effects that are often overlooked in the literature. 
First of all we introduce the effect of relative velocity of dark 
matter and ba ryonic fluids after recombinatio n. This effect, 
first studied in iTseliakhovich fc Hiratal (|2010l '). is nominally 
a second order effect in the perturbation theory of density 
evolution and hence has been ignored in studies based on 
the linear theory. Secondly, we emphasize the importance of 
a correct treatment of the sound speed variations in the time 
between recombination and z ~ 200 due to residual Comp- 
ton heating of the electrons on the CMB photons. As we 



show later in the paper, both effects play a significant role 
during the epoch of first halo formation and dramatically 
impact gas fractions in the first bound objects. 



2.1 Relative velocity of dark matter and baryonic 
fluids 

Before recombination, the baryons are tightly coupled to the 
photons and the sound speed is highly relativistic, ~ c/\/3. 
As the universe expands and cools the electrons recombine 
with the protons and the universe becomes transparent, 
leading to a kinematic decoupling of the baryons from the ra- 
diation around ^doc ~ 1000. After recombination the sound 
speed of the baryons drops precipitously down to thermal 
velocities, whereas the dark matter velocity remains unaf- 
fected, and after recombination the dark matter motion with 
respect to the baryons becomes supersonic. The relative ve- 
locity can be written as: 



Vbc(k) = -[6(b(k)-6ic(k)] 



(1) 



where k is a unit vector in the direction of k, and 9 = a~^V- 
V is the velocity divergence (we use comoving coordinates). 
The variance of this relative velocity is 



('^bc(x)) 



A?(fc) 



k 



I 



(2) 



where A^(fc) = 2.42 x 10~^ is the initial curvature pertur- 
bation variance per Infc. Integration of Eq. © at the time 
of recombinatiorQ (zrec = 1020) shows that the dark mat- 
ter moves relative to the baryons with root-mean square 
velocity ~ 30kms~^ corresponding to a Mach number of 
M = Vbc/cs ~ 5. This supersonic relative motion allows 
baryons to advect out of the dark matter potential wells 
and significantly suppresses the growth of structure at wave 
numbers higher than 



aH 



M 



40 Mpc" 



(3) 



were fcj is the Je ans wave number. 

As shown in ITseliakhovich fc Hiratal (|2010l 'l the relative 
velocity of the baryons and cold dark matter (CDM) is co- 
herent over scales of several comoving Mpc and the velocity 
in each coherence region is well-described by a 3-dimensional 
Gaussian probability distribution with variance 



0"vb( 



= <!vbc(x)0. 



(4) 



(Note that this is the total variance, i.e., including velocities 
in all 3 directions; the variance per axis is smaller by a factor 
of 3.) 

To see how the relative motion of baryons and dark mat- 
ter affect the formation of the first objects we need to solve 
a system of evolution equations that incorporate this effect. 



^ Technically the effe ctive redshift of kinematic decoupling 
llEisenstein fc HuiHooS) . since recombination is an extended pro- 
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The system of equations describing a high-fc perturbation 
mode in the presence of a background relative velocity is 

dt 



95b 
dt 



and 



^TT ^ — (S2c0c + iibOb) - i-nPb H 5-Ob- 

OT 2 



(5) 



The Wbc terms are nominally second order in perturba- 
tion theory, and hence one may wonder why they, rather 
than other second-order terms, are included. The reason is 
that the expansion parameter for these terms is not the den- 
sity perturbation S, but rather the ratio of the advection 
terms (e.g. v'^^'^ k5a/ a in the 5c equation) to the linear terms 
(e.g. dSc/dt ~ Sc/H). This ratio is 

(6) 



aH 



One can see that this expansion parameter increases as one 
goes to smaller scales and is of order unity at fc ~ fcvbc. 
Thus the iibc terms become nonperturbative at small scales 
k > fcvbc, and when treating these small scales one must 
keep these terms even if they are formally higher order in 
the perturbation theory. 



2.2 Complete heating model 

The system of equations of Eq. ([5]) assumes a spatially uni- 
form sound speed whi ch is a good first-order a pproximation. 
However, as shown in lNaoz fc Barkanal (|2005l ). it underesti- 
mates baryon density fluctuations by more than 30 percent 
at 2 = 100 and 10 percent at z = 20 for large wavenum- 
bers. A fully correct treatment of baryon density evolution 
requires analysis of the Compton heating from the CMB on 
the sound speed a nd fluctuations in the te mperature distri- 
bution. Following iNaoz fc Barkanal (|2005l '). we re-write the 
sound speed term of the last equation of Eq. ([5)l as 



-5b 



{Sb + St: 



(7) 



where 5t is the temperature perturbation which evolves as: 



^ 2d5b ^ a^e(t)^- 

dt 3 dt t~, 



(8) 

The second term on the right-hand side accounts for the 
Compton scattering of the CMB photons on the residual 
electrons from recombination. Here Xe{t) is the electron frac- 
tion relative to the total number density of gas particle^. 



[2.725 K]/a is the mean CMB temperature, and 



,-1 _ 8 CTT c 
^7 — ^ 



:p0^ = 8.55xl0-"yr- 

me 



(9) 



^ This is different from the recombination literature, which often 
takes Xe to be normalized to the number of hydrogen nuclei. At 
low redshifts these differ by 8 per cent due to the presence of 
helium. 



where ctt is the Thomson scattering cross section, p° is the 
photon energy density at z = 0, and T is the average tem- 
perature of the baryons, which can be calculated using the 
first law of thermodynamics: 



dT 
dt 



-2HT- 



T) a~ 



(10) 



Accounting for Compton heating of the residual electrons 
by the CMB photons is especially important on small scales 
(fc > 1 Mpc~^), which are also impacted by the relative 
motion effect. 



3 FIRST HALOS AND THEIR GAS CONTENT 

Both of the effects discussed above have a significant impact 
on the evolution of density perturbations on small scales and 
affect the formation of the first dark matter halos, as well as 
the subsequent accretion of the baryons and the formation of 
the first stars. We investigate the specific effects by studying 
the change in the characteristic mass scale that divides gas- 
rich and gas-poor halos produced by the relative velocity of 
the dark matter and baryonic fluids. 



3.1 Filtering mass 

In the ACDM universe, virialized dark matter halos form 
hierarchically on extremely small scales at very early times 
and start accreting baryons into their potential wells. If ha- 
los are heavy enough, accretion proceeds to the point where 
baryons start cooling through molecular line emission, con- 
densing into the first stars and galaxies. This accretion is 
counteracted by the bulk motion of baryons with respect 
to dark matter as well as by the thermal gas pressure. The 
combination of the two effects leads to the presence of the 
minimal halo mass scale at which baryons are still able to 
effectively accrete onto a halo. 

To study the effect of halo formation and baryonic ac- 
cretion it is convenient to divide space into a large number 
of patches of the size of the relative velocity coherence scale. 
In each patch with a given mean density and bulk velocity, 
we follow the evolution of density perturbations including 
the spatial variation of the baryo nic speed of sound due to 
Compton heating from the CMB (|Naoz fc Barkana|[2005l 'l. 

By evolving the system of equations © with the cor- 
rect sound speed term of Eq. ((7]) in each patch, we calcu- 
late the baryonic and dark matter power spectra. Their ra- 
tio is constant on large scales (small fc), and drops at high 
fc du e to the suppre s sion of growth by the baryonic pres- 
sure. ICnedin fc HuH | 19981 ) originally defined a "filtering" 
scale (essentially a time-averaged Jeans scale) that they used 
to identify the largest scale on which the baryon fiuctu- 
ations are substantially suppressed compared to those of 
the d ark matter. We use the generalized definition from 
iNaoz fc Barkaiial (|20o3), in which the baryon-to-total ratio 
is expanded to linear order in k'^, and written in the follow- 
ing form: 



Stot 



— = 1 - p- +rLSS, 

ft, rp 



(11) 



where the total density perturbation Stot = fhSh + fdmSdn 
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Figure 1. Perturbation ratio (5(,/5tot vs. comoving wavenumber 
k evaluated at 2: = 20 for tlie cases of v^,^ = (solid curve), 
Vhc = Ifubc (dashed curve), and i>bc = 'i'^vbc (dotted curve). In 
all cases overdensities are isotropically averaged over the direction 
of k with respect to v^c. 



(in terms of the cosmic baryon and dark matter mass frac- 
tions /b and /dm), and the fc-independent tlss term (which 
is negative) describes the relative baryon-to-total difference 
in the limit of large scale structure, i.e., where both the 
Vbc effect and the thermal pressure of the gas are negligi- 
ble (and restricted also to scales below the baryon acoustic 
oscillations) . 

In Fig. [T] we plot the isotropically averaged perturbation 
ratio Sb/Stot by averaging over the direction of k with respect 
to Vbc. On large scales the ratio is very close to constant, 
and using Eq. we can deduce tlss = —0.054 at z = 
20. The filtering scale kp is obtained by fitting Eq. (jlip to 
the calculated values of the ratio Sb/5tot from Fig. [T] This 
allows us to define the filtering mass in terms of the filtering 
wavenumber: 



47r / TT 



(12) 



where po is the mean matter density today. We n ote that 
this re lation is i of the definition originally used bv lGnedlnl 
l|200(]D . who also used a non-standard definition of the Jeans 
mass. 

The filtering mass plays an extremely important role in 
understanding the evolution of the first halos, as it provides 
a good approximation for the boundary between the gas-rich 
halos and halos that do not contain substantial quantities 
of gas. Traditionally one would assume that the separation 
between gas-rich and gas-poor objects is represented by the 
Jeans scale, which is the minimum scale on which a small gas 
perturbations will grow due to gravity overcoming the pres- 
sure gradient. However, the Jeans scale is related only to the 
evolution of the perturbations at a given point in time and 
does not account for significant variation of the Jeans mass 
with time. The filtering mass on the other hand reflects the 
baryonic pressure effects integrated over the entire history 
of the Universe, and provides a much better approximation 
to the boundary between gas-rich and gas-poor halos. 

An extensive study of the filtering mass properties 
and evolution history in models without the relative ve- 
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Figure 2. Evolution of the filtering mass with redshift in the 
regions with iibc = (thin solid line), tibc = Ifutc (dashed line), 
i;bc = 2(T„jc (dotted line) and global average case (thick solid 
line). We also show the evolution of the Jeans mass Mj (dot- 
dashed line). 



locity effect was p erform ed in iNaoz fc Barkanal l|2005h and 
iNaoz fc Barkanal l|2007l ). The properties of the filtering 
mass, however, are significantly modified in the regions 
where the bulk motion of baryons with respect to dark mat- 
ter potential wells is significant. In regions with high values 
of Vhc baryons tend to advect out of the collapsing dark 
matter halos, significantly increasing the filtering mass. We 
demonstrate this in Fig. [5] where we plot the evolution of the 
filtering mass with redshift in the regions with Hbc/cvbc ~ 0, 
1, and 2. We also show the globally averaged case by inte- 
grating the filtering mass over the full probability distribu- 
tion of the relative velocity, given by|f| 



Pvbc{v) = 



2 

vbc 



3/2 



47rii exp 



vbc 



(13) 



As noted earlier, the variance per axis is cr^bc/S. 

In Fig. [2] we also compare the filtering mass with the 
Jeans mass defined as: 



47r 

Mj = — po 



(14) 



where kj = •\/2/3aJ//cs is the Jeans scale (defined by set- 
ting the right-hand side of Eq. ((5| to zero, without the rel- 
ative velocity term, and neglecting here the correction of 
Eq. (O). Fig. [2] shows that the filtering mass reaches a max- 
imum value at redshift z ~ 40 (and generally varies only 
slightly throughout the plotted redshift range), whereas the 
Jeans mass continuously decreases with time due to the drop 
in the sound speed of the gas as the Universe cools. 

The filtering mass represents a time averaged Jeans 
mass and hence it decreases at the low redshifts, however, 
right after recombination baryonic perturbations on small 



^ This is the distribution of the magnitude of v^c, where the vec- 
tor Vbc is the result of linear perturbations and hence is drawn 
from a multivariate Gaussian. It thus happens to be the same as 
the Maxwell-Boltzmann distribution, even though the bulk ve- 
locities of baryons have nothing to do with thermal motions of 
particles. 
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Figure 3. This figure sliows the filtering mass Mp as a function 
of the relative velocity of the dark matter and baryonic fluids at 
2 = 20 (solid line) and z = 40 (dashed line). 
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Figure 4. The lower panel shows the dependence of the filter- 
ing mass on the angle 9 between the directions of v^c and the 
wavevector k (solid line), and the isotropically averaged value of 
Mp (dashed horizontal line), in regions where titc = 2(T„ijc a^t 
z = 20. 



scales are highly suppressed and they only catch up gradu- 
ally, causing the filtering mass to increase from low initial 
values. We emphasize that in the regions with a high value of 
the relative velocity the filtering mass is significantly larger 
than in the regions with small values of Ubc and hence the 
formation of gas-rich objects in those regions proceeds quite 
differently than in the regions with Ubc ~ 0. The filtering 
scale and mass (from Eqs. and (|12|) ) in regions with 

varying values of Vhc are given in Table [1] and the filter- 
ing masses at z = 20 and z = 40 are plotted in Fig. [3] In 
Fig. |4j we also show the dependence of the filtering mass 
on the angle 9 between the direction of Vbc and that of the 
wavevector k in regions where v\^c = ZctjjJc at z = 20; the 
plot shows that the main contribution to the filtering mass 
comes from the regions where the wavenumber k and the 
relative velocity vector Vbc are parallel (or anti-parallel). 



3.2 Gas content of the first galcLxies and minihalos 

We now investigate the amount of gas that falls into early 
haloes, and how much of this gas is capable of cooling. Here 
we use analytical approximations - namely the relation be- 





P{> ''be) 


kp (Mpc-i) 


Mp (Mq) 


4 


2.1 X 10""^ 


85 


7.75 X 10^ 


3 


5.9 X lO"*^ 


113 


3.37 X 10« 


2 


7.4 X 10-3 


166 


1.07 X 10*^ 


1 


0.392 


298 


1.85 X lO'^ 





1 


591 


2.39 X Iffi 



Table 1. Filtering scale and filtering mass for the isotropic aver- 
aging of the direction of k with respect to v^c a.t z = 20. 

tween the gas mass fraction /g and filtering mass Mp, and 
the Shcth & Tormen ( 1999) mass function - that have been 
calibrated against simulations with statistically isotropic ini- 
tial conditions and no bulk relative velocity. In our case with 
'i'bc 7^ the power spectra are both reduced and slightly 
anisotropic, but we expect these approximations to still be 
a useful guide since statistical anisotropy (e.g. 6'-dependent 
filtering mass) can only appear at second order in scalars 
such as the halo mass function or gas content. 

There is no a priori reason to suppose that the filtering 
mass, which is defined based on linear perturbations, can 
also accurately describe properties of highly nonlinear, viri- 
alized objects. Qualitatively we may argue that if pressure 
significantly opposes gravity during the halo formation pro- 
cess (which for some time is accurately described by linear 
theory) then it will significa ntly suppress t he amount of gas 
in the final virialized halo. iGnedinI (|200G| ) suggested based 
on simulations during cosmic reionization that the filtering 
mass accurately fits the mass for which a halo contains half 
the mean cosmic baryon fraction /(,, and fitted the simula- 
tion results to the following formula: 

-3/a 



fg — fh,0 



,a/3 



(15) 



where /b.o is the gas fraction in the high-mass limit. 
In this function, a higher a causes a sharper transi- 
tion between the high-mass (constant /g) limit and the 
low-mass limit (assumed to be /g oc A/^). This for- 
mula has subsequently been found to agree with hy- 
drodynamic simulations jNaoz. Barkana fc Mesingeij l2009l : 



iNaoz. Yoshida. fc Barkanal l20ld)" if we set 



/b,o = fh(l + 3.2n,RR) (jBarkana fc Loebl I2OI0I ). and use 
the filte r ing m ass as defined in Eq. (|12|l (which differs from 
iGnedinI (|200G| ). as noted earlier). Thus, at each redshift, in 
each patch of the Universe we may calculate the local value 
of Mp and from it the gas fraction in halos of various total 
mass. In Fig. O we plot the gas fraction as a function of 
halo mass in regions with varying values of relative velocity 
at 2 = 20. It is clear that halos that would be gas-rich in the 
Universe with no Wbc effect become gas-poor in the regions 
where the relative velocity is high. We also see that on av- 
erage, introduction of the v^c effect significantly lowers the 
gas fraction in all halos with Mh < IO'AIq. 

In order to find the total amount of gas in galaxies, we 
must integrate over the halo mass function in each patch. 
We start with the Wbc = case. Standard models for halo 
formation are based on spherical collapse calculations, in 
which the linear overdensity must reach a critical threshold 
5c{z) for the corresponding region to form a collapsed halo 
at redshift z. The halo abundance depends on the statistics 
of fluctuations on various scales, which can be parameterized 



© 0000 RAS, MNRAS 000, 000-000 



6 Tseliakhovich, Barkana & Hirata 




10- 



10° 10^ 



10° 



10" 



Figure 5. Change in tlie gas fraction fg at z = 20 as a func- 
tion of halo mass for regions with t)bc = (thin solid line), 
*'bc = l^ubc (dashed line), v^,^ = 2(T^i,c (dotted line), jj^c = ^<^vbc 
(dash-dotted line) and isotropically averaged case (thick solid line, 
nearly coincident with the dashed line). 



by the function S{R), the variance of fluctuations in spheres 
of radius R {S is also a function of redshift). In general 
we define f{5c{z),S)dS to be the mass fraction contained 
at z within halos with mass in the range corresponding to 
S to 5 -I- dS. We convert between halo mass AI and (initial 
comoving) radius R using the cosmic mean density. The halo 
abundance is then 



dn 

1m 



M 



dS 
dM 



(16) 



where dn is the comoving number density of halos with 
masses in the range M to M + dM. In the model of 



iPress fc Schechtei] (|l974 ). 



2/2 



2^S 



(17) 



where f = 5c{z)/\fS is the number of standard deviations 
that the critical collapse overdensity represents on the mass 
scale M corresponding to the variance S. 

However, the Press-Schechter mass function fits numer- 
ical simulations only qualitatively, and in particular it sub- 
stantially underestimates the abundance of the rare halos 
th at host galaxies a t high redshift. The halo mass function 
of lSheth fc TormenI (| 19991 ), which fits numerical simulations 
much more accurately, is given by 



1 



(a'i/2 



-a'i/2/2 



(18) 



with best-fit parameters (|Sheth fc TormenI |2003 ) a' = 0.75 
and q = 0.3, and where the normalization to unity is en- 
sured by taking A' = 0.322. 

The critical density of collapse Sc{z) is independent of 
mass and equals 1.69 in the Einstein de-Sitter limit, valid 
over a wide range of redshifts. Its value decreases at low 
redshift due to the cosmological constant, but more rele- 
vant for this paper is the decrease at very high redshift due 
to the effects of the baryons and radiation. The decrease 
is by ~ 0.05(1 -|- z) per cen t from the Einstein de-Sitter 
value (|Naoz fc Barkanall2007l ) , a small effect that is however 
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Figure 6. Dependence of o"(fbc|Afh) on the relative velocity v\^^ 
at z = 20 for a fixed mass of the collapsed halos: = W"^ Mq 
(solid hne), = 10*' Mq (dot-dashed line), = lO^M© 

(dashed hne), and = 10'' (dotted line). 



greatly amplified by the fact that galaxies at the highest 
redshifts form in halos that correspond to very rare density 
fluctuations. 

In our approach we must calculate how the halo mass 
function varies in differe nt regions. The Press-Schechter 
model has been extended l|Bond et al.lll99ll ) to describe the 
variation of the halo abundance in regions of various den- 
sity, and we can generalize this to include the bulk veloc- 
ity by including the variation of the function S{R) with 
Vhc- To demonstrate this dependence we plot in Fig. [6] 

= a{vhc\Mh) as a function of the relative velocity at 
a fixed halo mass Mh- We see that the variance of the den- 
sity perturbations decreases with increasing Vhc, leading to 
a delay in the collapse of dark matter halos. We also find 
that the change in a is larger for halos of low mass since 
the power spectrum on scales much larger than the filtering 
mass is unaffected by the relative velocity. 

We can summarize the important effects of the change 
in halo abundance and halo gas content at a given mass 
scale by calculating various gas fractions at each redshift. 
This can be done by using the Sheth-Tormen mass function 
and accounting for the changes due to the relative velocity 
effect. We calculate the fraction of the total matter density 
in halos above a certain mass scale by 



. . ^ f"° M dn 
Jm^ Po dM 



(19) 



and the fraction of the baryon density contained in those 
halos using Eq. (|15l) : 



/gas(> Mh) 



M dn fg 
M, po dM ft ' 



dM 



(20) 



We plot both fractions in Fig. [7] at z = 20 for Ubc/o"vbc = 0, 
1, and 2, and for the globally averaged case, where we take 
into consideration the global distribution of the Wbc- The 
plot clearly shows that in regions with high relative velocity 
the gas fraction in halos is dramatically suppressed. The 
global average (which comes out very close to the ^bc/cvbc = 
1 case) gives a suppression by a factor of 2.7 of the total 
gas fraction in halos. In order to separate out the various 
effects, we plot one case in which we use the correct halo 



© 0000 RAS, MNRAS 000, 000-000 



Suppression and Spatial Variation of Early Galaxies and Minihalos 7 




Figure 7. The mass fraction in halos above (upper panel) 
and the gas fraction in halos (lower panel) at redshift z = 20 for 
the case iibc/i^vbc = (thin solid line), 1 (dashed line), 2 (dotted 
line), and for the globally averaged case (thick solid line). In the 
lower panel we also show the case where we fix the value of Mp 
as calculated for v^,^ = and use the correct globally averaged 
halo mass function including the variation with ii^c (dot-dashed 
hne). 

mass function (as it varies with iibc) but fix the filtering mass 
to the Vhc = value. We find that the suppression arises from 
comparable contributions from the change in halo numbers 
(about a factor of 1.8) and from the reduction in the internal 
halo gas fractions (about a factor of 1.5). 

Stars are understood to form at high redshift out of gas 
that cooled and subsequently condensed to high densities in 
the cores of dark matter halos. Since metals are absent in the 
pre-stellar universe, the earliest available coolant is molec- 
ular hydrogen (H2), and thus the minimum halo mass that 
can form a star is set by requiring the infalling gas to reach a 
temperature of several hundred Kelvin required for ex citing 
H2 to the J ^ 2 rotational levels (jTegmark et al.lll997l '). This 
has been confirmed with high-resolution numerical simula- 
tions containing gravity, hyd rodynamics, and chemical pro- 
cesses in t he primordial gas ( Abel, B ryan fc No rman 2002!; 
Bromm. Coppie fc Larson 2002; Fuller fc Couchman 200Qi; 
Yoshida et al.l l2003l : iReed et alj I2OO5I ). These simulations 
imply a minimum halo circular velocity Vc ~ 4.5 km s~^ for 
forming a star, where Vc = \J GM/R in terms of the halo 



virial radius R. The simulations show that in a halo above 
the minimum mass (which at z = 20 is Mmin « 6 x 10^ Mq), 
the gas cools in the dense center and forms at least one star 
very quickly; this is understood theoretically since both the 
cooling time and the dynamical time are a small fraction of 
the cosmic age at that time. We are thus interested in the 
total gas fraction in halos above this cooling threshold; if 
there is a fixed star formation efficiency in these halos, then 
the this gas fraction is directly proportional to the stellar 
density in each region. 

We plot in Fig. [8] the evolution of various gas and to- 
tal mass fractions as a function of redshift. Even without 
the relative velocity effect, there is some (spatially-uniform) 
suppression predicted for the gas fraction in halos that 
can cool (i.e., a suppression of the overall star formation) 
by a factor of 1.2 at z = 20 and 1.5 at z = 40 (rela- 
tive to the cosmic baryon fraction); this is due to the fact 
that the baryon perturbations are still catching up to the 
dark matter perturbations at these redshifts, even on large 
scales (beyond the filtering mass), and simulations suggest 
that non -linear halo formation amplifies the remaining dif- 
ferences (^ Naoz. Yoshida. fc Barkanall2O10l : iBarkana fc Loebl 
I2OIO). The relative velocity effect adds an additional sup- 
pression of cosmic star formation by a factor of 1.6 at z — 20 
and 3.4 a.t z — 40. The relative velocities have a larger effect 
on the gas in minihalos, the smaller halos that accrete gas 
that cannot cool. Since the total mass fraction in halos con- 
tinues to increase as we consider smaller and smaller halos 
masses (Fig. [7|, the total amount of gas in minihalos is very 
sensitive to the filtering mass, which is what produces the 
(gradual) low-mass cutoff in gas accretion onto halos. In the 
absence of the relative velocities, the total gas fraction in 
halos at 2; = 20 is 1.5 x 10~^, consisting of 1.0 x 10~^ in 
minihalos and 5 x 10^'^ in galaxies. At 2 = 40, these gas 
fractions are 2.4 x 10"^, 2.3 x 10~^, and 1 x 10"*^, respec- 
tively. The relative velocities, in the global average, reduce 
these fractions to 6 x 10"^, 3 x 10"^, and 3 x 10"^ at 2 = 20, 
and 1.5 x 10-^ 1.1 x 10-^ and 4 x lO"'' at 2 = 40. Note 
that the gas fraction above the H2 cooling mass is really an 
upper limit to the gas fraction that undergoes star forma- 
tion. Any significant feedback effect will raise the effective 
threshold for star formation, making the total gas fraction 
in halos correspond almost completely to star-less halos (see 
discussion in Sec. 3)). 

3.3 Probability distribution functions 

In addition to plotting mean values of various quantities, it 
is interesting to consider their variation in various regions. 
Above, we have explicitly varied Wbc but averaged over the 
density fluctuations (note that the density and velocity are 
uncorrelated at a given point). However, in order to calculate 
the full amount of variation of various quantities, i.e., the 
probability distribution function (PDF), we must explicitly 
vary both the mean density and the value of Vhc in each 
region. 

In the absence of relative velocities, the extended Press- 
Schechter model gives the variation of the Press-Schechter 
halo mass function in regions of various mean density. No 
analytical generalization of this formalism is known for the 
more accurate Sheth-Tormen model, but iBarkana fc Loebl 
(2004 ) suggested a hybrid prescription that adjusts the 
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Figure 8. In the left panel we plot the total mass fraction in halos above the cooling mass (dashed lines) and total gas fraction in halos 
above the minimum cooling mass (solid lines). In the right panel we show total gas fraction in halos (dashed lines) and gas fraction in 
minihalos, i.e., in halos below minimum cooling mass (solid lines). All plots show two cases: no (regular lines) and the correct case 
where the effect is taken into account and the isotropic averaging is performed (thick lines). 



abundance in various regions based on the extended Press- 
Schechter formula, and showed that it fits a broad range of 
simulation results. Generalizing this prescription to include 
the effect of relative velocity, we set 



fhi!,s{dc{z),5ii,R,M,vt,c) = 



f ST {Sc{z),S') 



fpsi5ciz),S')_ 



S'iR)) , 
(21) 



where the mean overdensity in the patch is Sr, and for a 
given halo mass, the variance S' is calculated using the 
power spectrum modified by the local bulk velocity. The 
subtraction of S'{R) accounts for the fact that Sr arises 
from density modes on scales larger than the patch size, 
leaving only the remaining variance S' ~ S'(R) for fiuctua- 
tion modes within the patch to supply the additional density 
needed to reach 5c{z) and thus form a halo. In our case, the 
patches in which we will compute the baryon collapse frac- 
tion PDF will be spheres of radius R = 3 Mpc (comoving). 
Note that if used to compute a mass function, the above for- 
mula gives the Lagrangian halo number density, while the 
Eulerian density is larger by a factor of 1 + Sr; however no 
such transformation is necessary to compute the local frac- 
tion of gas in halos. 

We start by calculating the PDF for the filtering mass 
Mp. In the scenario without Ubc we would have a universal 
value of Mf, however, since various regions of space have dif- 
ferent values of relative velocity of baryonic and dark matter 
fiuids this produces a variation in Alp- The distribution of 
relative velocities is given by Eq. (|13|) and it translates into 
the distribution of Mf using: 



^ cIMf' 



(22) 



The PDF of the filtering mass at z = 20 and z = 40 is 
plotted in Fig.l^l These distribution functions are essentially 
determined by the distribution of the relative velocity and 
exhibit clear peaks which correspond to values of the filtering 
mass around the maximum of the Wbc distribution, which 
occurs at Vhc ~ 1.2<jvbc- As noted before (Figs. [2] and [3]), the 
filtering mass does not vary much in this redshift range, but 
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Figure 9. The PDF of the filtering mass Mf at z = 20 (solid 
line) and z = 40 (dashed line). We consider cells of radius i? = 3 
Mpc, and include the variation of Vy^^ as well as the mean density 
in each cell. 



at z = 20 it is slightly more sharply peaked while the PDF 
at 2 = 40 extends more towards high values of M f ■ Filtering 
mass has a rather significant scatter with the full width at 
half maximum ~ 1.5 X IO'^Mq at the redshift oi z = 20. 

To better understand global properties of the first ob- 
jects we calculate probability distribution functions of the 
total gas fraction in halos as well as the gas fraction in ha- 
los above the minimum cooling mass Mmin- As noted earlier, 
these gas fractions are affected by the distribution of relative 
velocities as well as the distribution of large scale overdensi- 
ties 5r. We consider the PDF of gas fractions inside spher- 
ical regions ("cells") of radius i? = 3 Mpc, which are small 
enough that Wbc can be treated as roughly constant over a 
cell. We obtain the PDFs by running a Monte Carlo sim- 
ulation that generates random values of iibc using Eq. (|13|l 
and of the large-scale overdensity within the cell 5r using a 
Gaussian of variance S'{R). In Fig. [10] we show the PDFs 
of the gas fractions in halos above and below the minimum 
cooling mass Afmin « 6 x 10^ Mq at 2; = 20. We also show the 
same distributions for the case with no Vhc effect. The figure 
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Figure 10. The PDFs of the gas fraction in halos above the 
minimum cooUng mass (soUd line), and the gas fraction in mini- 
halos, i.e., halos below the minimum cooling mass (dashed hne) at 
2 = 20 for the case with the tJ^c effect (upper panel) and without 
the effect (lower panel). 



shows how minihalos would be dominant at z = 20 (by a 
factor of 2 compared to galaxies) , but since Ubc has a larger 
effect on the minihalos, it makes the gas content roughly 
equal between galaxies and minihalos at this redshift. Each 
PDF has a non-Gaussian extension towards high fractions 
(in fact, the distribution is approximately lognormal). Thus, 
the peak of the PDF is significantly lower than the mean 
value; without the relative velocity it is 0.002 for galaxies 
and 0.008 for minihalos, and Vhc moves it to ~ 0.0015 for 
both. Also, the relative velocities reduce the full width at 
half maximum from 0.004 (galaxies) and 0.008 (minihalos) 
to 0.003 for both. 



4 CONCLUSIONS 

We have shown that the relative velocity of baryons and 
dark matter has a significant impact on the properties of 
the first bound objects and has to be considered in detailed 
studies of the epoch of reionization and especially earlier 
epochs. The supersonic motion of the baryonic fiuid relative 
to the underlying potential wells created by the dark matter 
causes advection of small-scale perturbations by large-scale 
velocity fiows, leading to a significant suppression of gas ac- 
cretion during halo formation and dramatically increasing 
the characteristic mass of gas-rich objects at high redshifts 



(z > 10). In particular, instead of this characteristic filter- 
ing mass being close to the Jeans mass of 2 x 10* Mq at 
2 = 20, it varies among various regions from this value up 
to ~ 10® M0, with a la value (and global average) around 
Mf = 2 X IO^Mq, i.e., an order of magnitude higher than 
without the relative velocity effect. 

The relative velocity effect also modifies the star forma- 
tion history, delaying star formation and causing significant 
spatial fiuctuations. However, since the minimum mass for 
H2 cooling 6 X IO^Mq at z — 20) is somewhat higher 
than the average Mf, the suppression effect of «bc is lim- 
ited to about a factor of 1.6 at z = 20 (added on top of 
the spatially-uniform factor of 1.2 from the still-depressed 
baryon perturbations on large scales), compared to a much 
larger effect (a factor of 3.3) on the gas fraction in star- 
less gas minihalos. The importance of the relative velocity 
grows steadily with redshift, so that at z = 40 the suppres- 
sion factors due to Vhc increase to 2.5 for galaxies (on top of 
a pre-existing factor of 1.5) and 21 for minihalos. 

In our detailed treatment, we included the spatial varia- 
tion of the baryonic sound speed, the suppression of baryonic 
perturbations on large scales, and the effect of the relative 
velocity, through the modified power spectrum, both on the 
halo mass function and the internal gas fractions in halos. 
In order to gauge the induced spatial variability, we further 
calculated the full probability distribution functions of the 
characteristic mass and of gas fractions inside of the first col- 
lapsed halos. These results are important for understanding 
of the relative velocity effect on large scales, and we plan to 
study them further. 

Our results significantly e xtend the work done recently 
bv lDalal. Pen fc SeliakI (|2010l '). For example, we find a sup- 
pression factor due to v^^ at z = 20 of 1.6 and 3.3, for 
star-for ming halos and minihalos, respectively. In their ap- 
proach iDalal. Pen fc SeliakI (|2010l ) did not separate these 
two categories, and found a factor of 2.5 suppression in the 
collapsed fraction, which under their approximation can be 
interpreted as a suppression of star formation. In our work 
we removed this and man y other approximations used in 
iDalal. Pen fc SeliakI l|201Cl ). Comparing to our work, we ex- 
pect that their calculation of Lyman-a fiux fluctuations is 
qualitatively correct but may be somewhat overestimated 
and requires a more detailed analysis. 

As we were finishing this paper, two simula tion papers 
appeared on the effect of Vhc at high redshift (jMaio et al.l 
I2OIOI : IStacv et aLlbOlol ). While both found a smaU suppres- 
sion of star formation, their results appear at first glance 
to show a smaller suppression effect than we predict. This 
difference is not surprising if we note that these simulation 
papers focused on star- forming halos around z ~ 15, while 
the largest effects that we find occ ur for star- l ess mi niha- 
los at higher redshifts. At z > 20, iMaio et al.l l|2010f ) find 
tens of percents difference in the gas fractions, although the 
statistical errors are large. ^Stacy et al. (2010) find a delay 
in gas collapse by Aa/a = 0.14 for Wbc/o"vbc = 1. We also 
note that the choice of initial conditions should be carefully 
considered: standard initial condition codes do not properly 
treat the separate baryonic and dark matter perturbations 
or the gas temperature perturbations, leading to a filter- 
ing mass that is too high by a factor of ~ 2 at Ubc ~ 
aoz. Yoshida. fc Barkanall2010l ): as such they may under- 
estimate the effect of relative velocities. 
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The simulations bv lMaio et all l|20ld) and IStacv et alj 
clearly represent a very important step and will serve 
as a good foundation for simulations with larger boxes and 
improved initial conditions. Eventually we hope that sim- 
ulations including Ubc will advance to the point where im- 
proved fitting functions for the local halo mass function and 
gas mass fraction become available. 

We note that various feedback effects may reduce or 
mask some of the effect of the relative velocity. For galaxies, 
local feedback from star formation may effectively raise the 
minimum halo mass for star formation (except for the very 
first generation of stars) . The possibilities include supernova 
feedback as well as radiative feedback acting via photoheat- 
ing and photoevaporation or suppression of H2 formation, al- 
though "positive" feedback due to X-ray ioni zation enhanc- 
ing H2 for mation has also been suggested (|Haiman et al.l 
1 19961 . 1 19971 ). For minihalos, astrophysical heating, e.g., from 
an early X-ray background, may heat the gas and raise the 
filtering mass above the value due to Ubc- There are many 
unknowns, but these various effects could begin to be sig- 
nificant hy z ~ 20, and very likely by the time of signifi- 
cant cosmic reionization. Still, the relative velocity between 
baryons and dark matter is the main determinant of the gas 
content of halos at the highest redshifts. 
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